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^ ' Abstract 



We present a consistent effective theory that violates the null energy condition (NEC) without 
developing any instabilities or other pathological features. The model is the ghost condensate with 
the global shift symmetry softly broken by a potential. We show that this system can drive a 
cosmological expansion with H > 0. Demanding the absence of instabilities in this model requires 
H < H^. We then construct a general low-energy effective theory that describes scalar fluctuations 
about an arbitrary FRW background, and argue that the qualitative features found in our model 
are very general for stable systems that violate the NEC. Violating the NEC allows dramatically 
non-standard cosmological histories. To illustrate this, we construct an explicit model in which the 
expansion of our universe originates from an asymptotically flat state in the past, smoothing out 
the big-bang singularity within control of a low-energy effective theory. This gives an interesting 
alternative to standard inflation for solving the horizon problem. We also construct models in which 
the present acceleration has w < —1; a periodic ever-expanding universe; and a model with a smooth 
"bounce" connecting a contracting and expanding phase. 
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1 Introduction 



According to our present understanding of cosmology the history of the universe has been character- 
ized by an incessant slowing down of the Hubble rate H = a/ a. Even the "accelerated" expansion 
(a > 0) today and during inflation at best corresponds to H = constant, and therefore can only halt 
the decrease of H. Intuitively, the fact that H cannot increase is a manifestation of the fact that 
gravity is generally attractive, and in some sense never "too repulsive." More formally, H < can be 
viewed as a consequence of the null energy condition (NEC), which requires that for all null vectors 
the matter stress-energy tensor must satisfies Tf^i^n^n'^ > 0. In a Friedmann-Robertson- Walker 
(FRW) spacetime this condition reduces to p + p > 0, which for a spatially flat universe directly 
implies H <0. 

But if the NEC is satisfied and H is always decreasing, why is the universe expanding in the 
first place? The conventional view is that H increases going backward in time until we reach an 
era where H ~ Mp\ (the "big bang") where quantum gravity effects become important. At this 
scale general relativity breaks down as an effective theory, so in this view the mystery of the origin 
of the expansion of the universe is inseparable from the issue of the UV completion of gravity. On 
the other hand, if the NEC is violated it is possible that the history of the universe is drastically 
modified without quantum gravity effects becoming important. For example, the present expansion 
might be due to a smooth "bounce" from an earlier contracting phase. An even more interesting 
possibility, which is just a limiting case of a bouncing cosmology, is a universe that asymptotically 
in the past is Minkowski space and as time goes on it simply starts expanding connecting smoothly 
to the subsequent familiar FRW cosmology and thus resolving the Big Bang singularity. These 
possibilities are discussed more fully below. 

Another motivation to consider violations of the NEC come from observations of the current 
accelerated expansion of the universe. These are conventionally summarized by the equation of 
state parameter w = p/p. A cosmological constant has w = —1, but present data allow w > —1.2 
[Tl[2]. The NEC implies w > —1, and it is natural to ask whether w < —1 is a viable possibility. 

One simple way to violate the NEC is to add fields with wrong-sign kinetic terms. Scalar 
excitations have energy unbounded from below (they are "ghosts"), and are therefore subject to 
catastrophic instabilities at all scales. In particular the vacuum is unstable to decay into ghosts 
and gravitons with an infinite rate in any theory that is Lorentz invariant in the UV [31 [1] . This 
kind of model can nonetheless be made compatible with observation by coupling it only to gravity 
and postulating that unknown UV physics breaks Lorentz invariance and cuts off the divergence 
for short wavelenths [Sj Hj. In this approach, the stability of the universe depends on unknown UV 
physics as the leading instability is <;t .hort distances 

The NEC is believed to hold for all well-behaved systems without instabilities. Other energy 
conditions can be violated by a sensible system simply by adding a suitable (positive or negative) 
cosmological constant. The NEC cannot, since it is saturated by a cosmological constant. Indeed it 
has been shown in Ref. [6] for a very broad class of models that whenever the stress-energy tensor 
violates the NEC the system has catastrophic instabilities, either in the form of ghosts or in the 
form of exponentially growing modes with arbitrarily short wavelengths ("tachyons")H Both types 

""^See [S] for an explicit realization of a Lorentz violating cut-off of a ghost instability. 
^The only exceptions are certain anisotropic systems which are stable but admit super luminal excitations 
[B]. However such exceptions are not relevant for cosmology, where one is only interested in isotropic systems. 
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of instability occur for all wavelengths down to the UV cutoff of the theory. Therefore, unlike the 
Jeans instability, these instabilities cannot be damped by Hubble friction, which is only effective at 
frequencies smaller than H. 

The results of Ref. [6] are derived at the two-derivative level in a systematic low-energy expansion. 
Higher derivative corrections are suppressed by some scale M, and in general they can be safely 
ignored at frequencies and momenta smaller than M. However this is not the case if the two- 
derivative action is degenerate, and some modes propagate only thanks to higher derivative terms. 
This is precisely what happens in the ghost condensate [7] and in more general theories of massive 
gravity [8]. For example, scalar fluctuations around the ghost condensate have dispersion relation 
up' = k'^/M'^. Theories of modification of gravity can be thought as non-trivial scalar backgrounds 
with a stress-energy tensor equal to the one of a cosmological constant ^ . Therefore they lie on the 
threshold of violating the NEC. This explains the origin of the degenerate dispersion relations and 
tell us that these theories are a first step towards violating the NEC. 

In this paper we use the ghost condensate as a starting point to construct a consistent and 
technically natural low-energy effective field theory that violates the NEC without any instabilities. 
The simplest example is a ghost condensate (/> rolling up a linear potential, which gives rise to 
superacceleration (H > 0). The late-time solution for (p climbs up the potential indefinitely, making 
the universe expand with an indefinitely growing Hubble rate, H cx \/t. This solution is an attractor, 
in the sense that trajectories flow to a common late time behavior, like in conventional slow-roll 
inflation. An interesting feature of this model is that we cannot violate the NEC by an arbitrary 
amount, in the sense that we must have 

H^H^ (1) 

to avoid instabilities. This comes from an interplay between gradient and Jeans instabilities. As 
discussed in Ref. [6], violating the NEC gives rise to a dispersion relation of the form ~ —k'^, 
which has a gradient instability. This instability is cured at short distances by the higher-derivative 
term giving a dispersion relation schematically of the form lo'^ ~ —k"^ + k'^, which is unstable only 
for long times u) < ti-'grad- However, the k^ term increases the mixing with gravity, leading to a Jeans 
instability similar to that for ordinary matter, so the system is unstable for lo < <-ij jeans- Remarkably, 
we find the model-independent relation between the instabilities 

^'-'grad Jeans ~ 

H. (2) 

Hubble friction will damp all instabilities only if both Wgrad and wjeans are smaller than H, which 
leads to Eq. ([T]). We emphasize that if this inequality is violated, the instability does not extend 
to arbitrarily short distances and time scales, and is under control of the effective theory. These 
instabilities are therefore much milder than ghost or tachyon instabilities, and may in fact have 
interesting cosmological consequences. (For example, the conventional Jeans instability gives rise to 
structure formation in the universe.) 

In fact, these features are much more general than this particular model. We demonstrate this by 
performing a general analysis of scalar fluctuations about an arbitrary FRW cosmological history. 
The analysis can be viewed most simply as the general case of an expansion driven by a rolling 
scalar. However, it also applies to the fluctuations of the Goldstone mode of time translations in 
the case where the expansion is driven by a mixture of fluids and rolling scalar flelds, and so the 
analysis is quite general. We find that the qualitative features of the ghost condensate in a linear 
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potential are completely generic. In particular, whenever H is positive the system is unstable at 
the two-derivative level, and the same interplay between gradient and Jeans instabilities leads to 
Eq. ([2|), so we must have H < in order for instabilities to be completely absent. 

We then consider a number of cosmological applications of models based on ghost condensation 
that violate the NEC without instabilities. We first discuss the possibility that the present expansion 
of the universe is superaccelerating (H > 0), corresponding to w < —1. The model consists of the 
ghost condensate rolling up a potential. In such a model H can in principle increase smoothly 
forever, although at some point H becomes so large that the model exits the regime of validity of 
the effective field theory. To get a measurable violation of the NEC we must be close to saturate 
the inequality ([T|), suggesting that there may be long wavelength growing modes in the dark energy 
sector that just started evolving. These may provide an additional signal for this class of models. 

We then construct a model where the universe starts from Minkowski space in the asymptotic 
past. The model consists of a ghost condensate with <j) = t that rolls up a potential V ~ (j)""^ ■ For 
large negative t, this model gives H ~ so the universe goes from a zero curvature, zero energy 

state to higher curvatures and energies, and this can then be smoothly connected to a standard 
FRW expansion. 

Next, we consider a model which gives rise to a cyclic universe that is always expanding. While 
the scale factor a{t) steadily grows larger and larger, the Hubble parameter and the energy density 
are periodic functions of time. In this model, a ghost condensate travels along a periodic potential, 
giving rise to a phase of superacceleration followed by reheating, a radiation-dominated and then 
matter-dominated phase, followed again by superacceleration, and so on. In such a scenario, the 
present accelerated phase is the beginning of primordial inflation! 

Finally, we consider the possibility that the universe smoothly bounces from a contracting to an 
expanding phase. Such a possibility has been considered previously in the literature (for a review, 
see Ref. [9]), but in previous treatments the bounce was not under theoretical control. We construct 
a completely smooth bouncing solution with no instability. Note that since H ^ {) aX, the bounce, 
while H ^ 0, there is always an unstable mode violating Eq. ([T|). However, it is easy to arrange 
that instabilities do not have time to grow much during the bounce, so that this phase is under 
theoretical control without catastrophic instabilities. 

We conclude that it is possible to have physical systems that violate the NEC without instabilities 
or other pathologies, and that this opens up a large number of interesting new possibilities for 
cosmology. The illustrative examples we consider are intended only as toy models; we leave the 
construction of realistic models and their phenomenology to future work. 



2 Superacceleration without Instabilities 

The simplest example of a system that violates the null energy condition without developing in- 
stabilities can be obtained as a deformation of the ghost condensate [7] , as outlined in [TO] . Ghost 
condensation can be realized starting from a Lagrangian for a derivatively coupled scalacl 

£ = ^ M^P{X) , X = -gf^'d^cPducI) (3) 
^We use a nietric with (— , +, +, +) signature. 
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with M an arbitrary mass scale. (Note that we take (f) to have mass dimension —1 so that X is 
dimensionless.) The absence of non-derivative couplings is natural if the scalar has a global shift 
symmetry 

(A ^ (/) + A. (4) 

In an expanding universe, one would expect that the field is asymptotically driven to rest {(j) 0) 
by Hubble friction. However, it is easily checked that there is also a cosmological solution at a 
minimum of P, 

(t) = ct, P'(c2) = 0, (5) 

where the metric is either Minkowski or de Sitter space. This is because at the minima of P the 
stress energy tensor is the same as for a cosmological constant, although the field is evolving with 
time [7J. Without loss of generality we can take c = 1 by a redefinition of the field (j). We assume 
that the vacuum energy is positive, so the metric is de Sitter space. However, unlike a cosmological 
constant, the ghost condensate has physical scalar fluctuations defined as 

(t){t, x)=t + 7r{t, x) . (6) 

The action for tt can be obtained expanding the original Lagrangian Eq. ([3|). Up to now we have 
neglected terms with more than one derivative acting on (p, like (Ocp)'^. Although they are not 
relevant for the unperturbed solution, they are important for the vr dynamics since they give the 
leading spatial kinetic term. We therefore obtain an action of the form 



S = d x^f^^ 



(7) 



where we have included a cosmological constant term and chosen P"(l) = \^ which amounts to a 
redefinition of the mass scale M. The action is not manifestly Lorentz invariant, as expected since 
the unperturbed solution spontaneously breaks the Lorentz symmetry. In particular there is no 
(Vvr)^ spatial kinetic term, but the leading term is (V^vr)^. The most generic action for vr can be 
obtained using only symmetry arguments: vr, as shown in Eq. ([6]), non- linearly realizes the broken 
time diffeomorphisms. This approach will be developed further below. 
To obtain a model with if > 0, we introduce a linear potential 

V = V'(j)^ V' = constant . (8) 

It is natural to have |y | <C because the potential is the only term that breaks the shift symmetry. 
In fact, in the absence of gravity a linear potential does not break the shift symmetry at all, since 
C ^ C + constant, and so the only radiative corrections come from gravity. These are very small 
given the low cutoff of the effective theorj0, so this form of the effective theory is technically natural. 

In the presence of the potential the background solution changes slightly, and it can be described 
by a homogeneous mode '7ro(t) for the field vr. Its equation of motion is 



^Constraints on the cutoff scale, M < 100 GeV, come from limits on IR modification of gravity in ghost 
condensation [71 lllj. 
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Neglecting the tt term, we find the attractor solution 

v 



3M^H 



The vf term is negligible only if 



H < 



(10) 



The velocity of the background is slightly reduced by the tilt as the field rolls up the potential. For 
small tilt, the stress energy tensor remains close to the one of a cosmological constant with a slowly 
increasing magnitude, so that H grows with time. Notice that the ttq perturbation becomes smaller 
and smaller with time as H increases and we approach the following asymptotic behavior 



3M2j 



oc t 



v 



oc 



3H{t)M^ 



(12) 
(13) 



In this solution, H constantly grows with time, while ttq approaches the minimum of the function P 
(see Fig. [I]). 
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Unstable 
gradient 
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Stable 
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Figure 1: Schematic representation of the model presented in section [2j 



We now proceed to study the stability of the system. First of all, we notice that because we 
are dealing with an accelerating background where modes exit the Hubble horizon and freeze out, 
an unstable mode has time to grow only if its rate is much faster than H. We can thus restrict to 
consider wavelengths and timescales much smaller than . This means that H and H do not 
change during the time scales of interest, and can be treated as constants. 

The quadratic action for vr can be obtained expanding Eq. ([7]) around the new background ttq: 




where we have used that at late times Y' ~ 6i7iirAfpj|£| The background ttq induces a (VTr)^ term, 
and for > it has the "wrong" sign, corresponding to P' < 0. This signals the presence of 

^There are small shifts of the coefficients of tt^ and (V^tt)^ that can be safely neglected. 
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instabilities, which we wih refer to as 'gradient instabihties'. In the absence of a tilt the expansion 
of the universe drives the field velocity to the point P' = 0, where there is no (V7r)2 term. This 
point separates the stable region P' > — where the fluctuation vr has a positive gradient energy — 
from the unstable one P' < 0. The negative tilt drives (j) into the unstable regime (see Fig. [1]). 
Neglecting the mixing of the vr mode with gravity, we obtain the dispersion relation 

+ M4 ^^^^ 

Note that stability is ensured for sufficiently short wavelengths by the k'^ term. The modes with the 
fastest instability rate are the ones with /c^ ~ Mp^H / M and their rate is of order 

• (16) 

These instabilities are absent if this rate is slower than the Hubble expansion, which requires the 
following constraint on the model parameters: 

The mixing of the tt mode with gravity gives rise to a second kind of instability, which is already 
present in the ghost condensate model in the absence of a potential. It comes from the mixing of 
the scalar with gravity and it can be interpreted as a sort of Jeans instability [7\. The instability 
rate in Minkowski space is of order 

'^Jeans ~ " ( ^^7^ ) • (18) 




Also in this case the instability is cured if the Hubble rate is sufficiently fast 

MM2 



Comparing Eqs. (|16p and (jlSp . we see that 

^grad ^Jeans ~ 

H , (20) 

independently of the model parameters. The two conditions for stability push in opposite directions, 
and putting them together we are left with the window: 

We conclude that at least for H <C the model parameters can be chosen in such a way that 
there are no instabilities, while H grows indefinitely. Notice that if Eq. (I21j) is satisfied at the initial 
time, it remains valid forever as H grows, while H /H decreases with time. 
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3 Effective Theory for General FRW Models 



We now turn to a much more general analysis of cosmological models, and show that the features 
found in the previous example are quite general. In particular, H < is always required to avoid 
exponentially growing modes in the class of models we consider. 

We would like to know whether a given cosmological expansion a{t) suffers from instabilities. 
Clearly, we need to know all the light degrees of freedom to answer this question, and we cannot 
answer it in complete generality. We will consider a framework which we believe is very general. 
Most simply, it can be described as the effective theory for fluctuations about a FRW background 
driven by a single rolling scalar with a completely arbitrary Lagrangian. This therefore includes for 
example a standard slow-rolling inflaton and all possible deformations of the ghost condensate, such 
as the one considered in the previous section. 

In fact, the analysis is even more general. We are effectively focusing attention on a scalar 
excitation which is present in virtually all expanding universes independently of what matter is 
actually driving the expansion: the perturbation corresponding to a common, local shift in time for 
all the matter fields tpm- That is, given a background homogeneous FRW solution a{t), xpmit), we 
consider the perturbation 

6i>m{x) = i^m{t + TT{x)) -i^m{t) , (22) 

parametrized by vr(x), and the corresponding scalar perturbation of the metric as imposed by Ein- 
stein's equationslfl This is the perturbation which in the long-wavelength limit is called adiabatic and 
obeys a conservation law which is insensitive to the matter content of the universe [l^ • By its very 
definition such a perturbation can be gauged away from the matter sector by an x-dependent shift 
in time. In other words, as long as we are only interested in this specific perturbation we can always 
choose a gauge in which the matter fields are unperturbed, 5iprn = 0, and the scalar fiuctuation is in 
the metric. We will adopt this gauge choice and refer to it as 'unitary gauge'. (For the case where 
the expansion is driven by a single scalar field, this corresponds to the gauge choice (p{t,x) = (j)Q{t), 
where (j)o{t) is the unperturbed scalar solution.) 

The existence of such a fluctuation mode sounds like a trivial consequence of the presence of 
a time-dependent FRW background, but it is not. If the expansion of the universe is driven by a 
mixture of rolling scalar fields, then the fluctuation we are interested in is indeed deflned by Eq. (j22p : 
it is the Goldstone boson associated to the spontaneous breakdown of time-translations. However in 
general the situation is more subtle. For instance the ground state of a solid or a fluid is characterized 
by three scalar condensates {cj)^{x)) = x* (z = 1,2,3) that spontaneously break the product of 
spatial translations and internal shift symmetries down to the diagonal subgroup [6]; in particular 
time translations are unbroken. Still, in the case of a fluid it can be shown that non-vorticous 
excitations — i.e., sound waves — can be described at the classical level as Goldstone bosons of broken 
iime-translational invariance around a rolling scalar condensate {4>{x)) = t, like a ghost-condensate 
[6]. So for any mixture of cosmic fluids and rolling scalars there exists an excitation that behaves 
like the Goldstone of broken time-translations, whose action is invariant under spatial translations 
and rotations. When coupling to gravity is taken into account, this translates into invariance under 

^Of course one such perturbation is trivial, corresponding to the FRW background solution with unper- 
turbed matter fields expressed in an unconventional set of coordinates. We eliminate this pure gauge mode 
by fixing a gauge for the metric. 
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time-dependent spatial diffeomorphisms + t). For more generic systems, e.g. solids 

or Lorentz-breaking massive gravity models [8], the low-energy degrees of freedom and the residual 
symmetries will be different, and our conclusions will not necessarily apply. Notice moreover that 
the presence of additional modes besides the one we are considering could, by mixing with it, change 
its dynamics. 

We now proceed to construct the most general effective action for the Goldstone, around a 
generic FRW background a{t). As we discussed, in unitary gauge the scalar mode does not appear 
explicitly in the action but it is part of the metric. Once we choose unitary gauge in fact time 
reparameterizations are not allowed anymore so that the metric contains an additional scalar degree 
of freedom. As we will discuss later the vr dependence of the action can be restored using the usual 
Stiickelberg procedure, as in the case of massive gravity [13]. In unitary gauge the full dynamics is 
described by an action for gravity which does not have the full diffeomorphism invariance, but it is 
only invariant under time-dependent spatial diffeomorphisms ^ + ^*(x, t). 

It is particularly convenient to work with ADM variables [H]. These are the 'lapse' N = 
the 'shift' Ni = goi, and the induced metric gij on hypersurfaces of constant t. In the 
following we will lower and raise spatial indices with the three-dimensional metric gij and its inverse 
g'^^ . In ADM variables the full 4D metric reads 

ds"^ = -N^dt^ + gij{dx' + N'dt){dx^ + NHt) . (23) 

The ADM formalism keeps manifest the invariance under 3D space diffeomorphisms: only quantities 
which are manifestly covariant under these transformations appear in the equations. The invariance 
under time reparameterization, although not manifest, is obviously still there. For our system, in 
unitary gauge the time variable is set by the unperturbed matter fields tpmit) so that the splitting 
between time and space, which is in general arbitrary, takes here a physical meaning. Therefore the 
unitary gauge action is not invariant under time diffeomorphisms, while all the unbroken symmetries 
(time-dependent spatial diffeomorphisms) are manifest in the ADM language. 
The Einstein-Hilbert action is expressed in the ADM language as 

Ser = ^M^i ld^x^R= ^Mli I d^x dt [iVi?(3) ^ l-(E'^Eij - E\^)] , (24) 

where R^^^ is the Ricci scalar of the induced 3D metric. Eij is related to the extrinsic curvature Kij 
of hypersurfaces of constant t, 

E,j = NK,j = ^ [dtOij - ^iNj - V,iV,] , (25) 

where V is the covariant derivative associated to the induced 3D metric gij. 

The full action consists of the Einstein action plus matter terms. The Einstein action is invariant 
under time reparameterizations, but the matter action is not because we have chosen the gauge vr = 
where the scalar fluctuations are parameterized by the metric. We expand about a given cosmological 
background 

ds^ = -dt^ + a^{t)dx^ , Vm = V'mo(i) (26) 
and consider general metric fluctuations 

N = l + 5N, Nj=6Nj, gij = a\t)6ij+5gij . (27) 
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The action written in terms of the fluctuation fields will have time-dependent coefficients because 
they are functions of the background matter fields ipmoit)- We therefore take the matter action to 
consist of the most general Lagrangian invariant under spatial diffeomorphisms, with time-dependent 
coefficients. Although the natural integration measure is \/—g d?x dt, we find it more convenient to 
use the 4D invariant measure \/—g d'^x = ^/gNd^xdt. Since is a scalar under spatial diffeomor- 
phisms, this simply amounts to a reshuffling of the terms in the Lagrangian. 

3.1 Two-Derivative Goldstone Action 

In this subsection, we use the formalism described above to construct the two-derivative action for 
the vr mode described in the previous subsection. At short distances and times where the mixing 
with gravity can be neglected, this describes the physical scalar fluctuation of the system. The 
results obtained here will be confirmed by a complete analysis of the full gravitational action in the 
section [3^ but the present analysis is much simpler and more transparent. 

Our strategy is to write the most general matter action in unitary gauge vr = order by order 
in the metric fluctuations 5N , 6Nj, and 6gij. We can then restore the vr dependence using the 
Stiickelberg trick. That is, given a space-diff invariant term in unitary gauge, e.g. l/N"^, we write 
its transformation law under time-diffeomorphisms 

^,=r^j^,+2d,e-{d,e)\ (28) 

where we kept terms of zeroth order in the metric fluctuations and of second order in the transfor- 
mation parameter . We can then restore invariance under time reparameterizations by introducing 
the field vr transforming as 

TT TT - ^° , (29) 

and using the gauge-invariant combination 1/A^^ -|- 2 5ovr — (S^vr)^ in place of 1/A'^^ wherever the 
latter appears in the unitary-gauge Lagrangian. This procedure, where we neglected the gravitational 
perturbations, gives the action for the Goldstone vr. In order to expand the action at quadratic order 
in vr we need the transformation law of all terms in the action under time diffeomorphisms at second 
order in the transformation parameter . The Einstein-Hilbert action is invariant under a generic 
4D diffeomorphism, thus from there we get no contribution to the vr Lagrangian. The same holds 
true for the measure \/—g d'^x in the matter action. Also all coefficients that explicitly depend on 
time should be evaluated at {t + n). However it is reasonable to assume that the typical time scale 
for these coefficients will be of order so that, after expanding in tt, these give rise only to non- 
derivative terms suppressed by H, H, etc.: we can thus neglect these contributions as long as we are 
interested in fluctuations with wavelengths much shorter than the cosmological horizon. Finally, we 
are interested in constructing the action for vr in a systematic derivative expansion. In unitary gauge 
this corresponds to an expansion in derivatives of the metric fluctuations, which we now proceed to 
write down. 

In order for the given FRW background to solve Einstein's equations, the matter action must 
contain 'tadpole' terms, i.e. terms that start linear in the metric fluctuations. Note that the linear 
terms are canceled by a tadpole term in the Einstein action, since we are expanding about a solution 
of the full action. However, the Einstein action is completely invariant under time reparameteri- 
zations while the matter action is not, so it is useful to keep the linear terms in the matter action 
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for the Stiickelberg trick. At the zero-derivative, hnear level the only invariant under spatial diffs 
is 6N. Apart from that, a linear piece in the fluctuation can also come from the \/—g that makes 
the integration measure diff invariant. Therefore we have two independent operators at this order. 
Instead of choosing the operators 5N and 1 as our 'basis', we find it more convenient to choose 
a different combination, and 1; the two bases are equivalent at linear order, their difference 

being quadratic in 6N; the reason of this choice will soon become clear. At this order the matter 
action is thus of the form 



'matter 



-9 



Ait) 



(30) 



As discussed above, we allow for generic functions of time c(t), A(t) as coefficients of the operators 
in the Lagrangian. In fact, the above action, truncated exactly at this level, is the complete action 
in the case of an ordinary scalar field (p with a potential V{(j)) in unitary gauge. To see this, just 
re-write the standard action for on a solution (f){t) using ADM variables: 



d^x 



-9 



1 



(90)2 - y(0) 



d'^x 



-9 



1 



2iV2 



0' - V{<Pit)) 



(31) 



which is exactly of the above form with c{t) = and A{t) = V{(j){t)). This is the reason why the 



parameterization in Eq. (j30p is convenient. 

In principle one could also add tadpole terms that involve derivatives of the metric, like for 
instance an extrinsic curvature term However by integration by parts such terms can always be 
rewritten as a combination of our non-derivative tadpoles plus derivative terms that involve higher 
powers of metric fiuctuations, which we will discuss later. 

We now rewrite the coefficients c{t) and A(t) in terms of the Hubble parameter and its deriva- 
tives using the Friedmann equations. The matter stress-energy tensor evaluated on the background 
configuration can be read from the action Eq. ()30p : 



2 SS„ 



-9 ^9 



Too = c{t) + A{t) 

Tij = a''{t)6ij{c{t)-A{t)) 



(32) 



Therefore Friedmann equations are 



3M2j 



a 

a 



3M2j 



[c{t)+A{t)] 
[2c{t) - A{t)] 



(33) 
(34) 



Note that quadratic and higher order terms in the matter action do not contribute to the background 
Einstein's equation, since their contribution to the stress-energy tensor is at best linear in the 
fiuctuations and therefore vanishes if computed on the background. Thus the above Friedmann 
equations uniquely determine the coefficients c{t), A{t). Eq. (j30p can therefore be written 



'matter 



d^X 



-9 



(35) 



We now make use of the Stiickelberg trick. Given the transformation law Eq. (j28p . we obtain simply 



matter 



d^X 



-g{Ml,H){d^) 



(36) 
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where we neglected the total derivative linear term 9o^lj 

The Lagrangian for it is that of a standard relativistic scalar, but with an overall —H factor. 
Note that whenever H > the kinetic term has the wrong sign, leading to a catastrophic ghost 
instability. 

We can try to cure the instability we just found by adding quadratic terms to the effective action 
Eq. ([35|l . At the zero-derivative level in ADM variables the only possible quadratic operator is 

^matter ^ ^matter + jd^X^ \M\t) {6Nf . (37) 

The coefficient M^(t) is unconstrained, and it can be a generic function of time. In realistic situations 
we expect its typical time-variation rate to be of order so we can approximate it as a constant 
when studying short- wavelength fluctuations. In order to reintroduce the Goldstone field we just 
need the trasformation law of at first order \n , N ^ N — Oq^,^. This changes the Goldstone 
action by 

As already stressed in Ref. |15j, for large enough — larger than Mp^H — the vr^ gets the healthy 
sign, but the gradient part of the action still has the wrong sign for H > 0, signaling the presence 
of exponentially growing gradient instabilities at all wavelengths. 

At the zero-derivative level in the metric, which corresponds to the two-derivative level in the 
Goldstone, there is nothing more we can do. Indeed in Ref. [6] it was shown in a broad class of 
theories that whenever the matter stress-energy tensor violates the null energy condition the system 
has instabilities, either in the form of ghosts or in the form of exponentially growing modes. This was 
shown at the two derivative level. The example of the last section makes no exception, its stability 
crucially relying on higher derivative terms in the Goldstone Lagrangian and on the presence of the 
horizon for cutting off the instability in the IR. We are therefore led to consider higher-derivative 
terms in the effective action. 



3.2 Higher-derivative Goldstone Action 

At next order in the derivative expansion we must include terms that involve the extrinsic curvature 
Kij of constant-t hypersurfaces. Of course, since N is invariant under spatial diffs, we can as well 
use its 'reduced' version Eij of Eq. (j25p . Notice that the background solution has non-zero extrinsic 
curvature, E^j^ = a^H gij, so we are interested in the fluctuation SEij = E^j — E^\ With it we can 
construct the quadratic operators 6E'^i6N, 6E^^ SEij, and SE^i^. Let us ignore for the moment the 
first of these operators, and let us concentrate on the other two, for which the discussion is simpler. 
In order to see what these operators correspond to in terms of the vr field, it is sufficient to consider 
the transformation law of 6Eij at linear order in the time-diff parameter and at zeroth order in 
the metric fiuctuation. From the trasformation law g^i, — > g^i, + V + ^ u^ti we get 

5Eij ^ 5Eij - 5ij dt{a^Hio) - did^Co ■ (39) 

'^This term is a total derivative as long as we consider momenta much larger than H. When the time- 
dependence of the coefficients is completely taken into account all linear terms are going to cancel exactly, 
since the background solves Einstein's equations. 
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The term proportional to 6ij is negligible for momenta much larger than H. In such a limit we have 



-9 



2' 



Jd^x 



-9 



(40) 



where M = + M'^, and as before we approximated the coefficients and M'^ as constant in 
time. 

According to our effective field theory approach, we are assuming that the scales M, M, and 
M' are all of the same order of magnitude, and all further terms appearing in the effective action 
will be weighted by the same scale. Then terms that involve higher time-derivatives on vr will be 
negligible with respect to Eq. (j38p at frequencies below this cutoff scale M. Similarly, terms involving 
more spatial derivatives on vr can be neglected with respect to Eq. (|40p at momenta smaller than 
M. Therefore in this regime the only possible other term in unitary gauge is a linear term in the 
3D curvature R^^\ However, as clear from Eq. (j24p such a term can always be re-written as an 
Einstein-Hilbert term plus terms quadratic in Eij. 

Collecting all of the terms considered so far, the matter action is 



I 



d^X 



fM* 



(41) 



The positivity of the vr^ term is guaranteed for — 2Mp^H > 0. For simplicity in the rest of this 
section we assume — 2Mp^H ~ M^. Then the dispersion relation of vr excitations is exactly the 
one we studied in the example of section [2] 



M4 



M4 



(42) 



where k is the physical momentum. For positive H, the system has oscillatory excitations with 
dipersion law oc at large momenta and exponentially growing instabilites at low momenta, 
cx — fc^. That is, the higher derivative term cures the instability at small wavelengths. The critical 
momentum below which the system is unstable is 



k, 



'grad 



M 



(43) 



corresponding to a typical instability rate 



PI 



MM2 



(44) 



We are interested in pushing the instability time-scale outside the cosmological horizon, where the 
Hubble friction freezes the dynamics of perturbations. Note that the critical length-scale can even 
be much smaller than the Hubble horizon — what matters is the typical time-scale for developing the 
instability. 

Up to now there is no obstacle in making the instability rate as slow as we like. It is enough to 

2 

choose a very large M coefficient for the higher derivative term to make tijgrad niuch smaller than 
H. But here is where the mixing with gravity becomes crucial. Usually, if we restrict to momenta 
and frequencies much larger than H we are allowed to neglect the fact that matter fluctuations are 
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mixed with gravitational ones. This is because, by definition of T^^u, such a mixing comes from a 
Lagrangian term of the form h^u 6T^'^ , where JT^'^ is the fluctuation in the matter stress-energy 
tensor at first order in the matter fluctuations. Usually this is non-zero only if the background 
stress-energy tensor is non-zero. For instance, for an ordinary massless scalar (j) one has 

h^^ST^"" r^h^odifr. Mpi Hhdcp , (45) 

where 0o(O is the background field configuration and (p is the fluctuation, and we assumed that the 
expansion of the universe is dominated by (p. The mixing term is suppressed by H, and at momenta 
and frequencies much larger than H it can be safely neglected. In the case of a standard fluid one 
can repeat the same estimate, for instance by means of the effective Lagrangian presented in Ref. [6]. 
The result is exactly the same: the mixing between sound waves and gravity is proportional to H. 
As a consequence, the Jeans-instability time for a fluid that drives the expansion of the universe is 
always of order H~^. In our case the situation is very different: we already stressed that all quadratic 
terms we add to Eq. (j35|) gives contributions to the matter stress-energy tensor that are linear in 
the fluctuations, weighted by coefficients that are unrelated to the background energy-density. This 
means that the mixing they induce is completely unrelated to H, and can lead to modifications of 
the fluctuation dispersion law at frequencies parametrically larger than H. 

In particular, the higher-derivative terms we are adding to cure the gradient instability enhances 

2 

the mixing of vr with gravity. The larger we make M in order to cure the gradient instability, the 
bigger the effect of this mixing, and — if the mixing tends to destabilize the system — the faster the 
gravity-induced instability. We will see in the next section that indeed when the full gravitational 
dynamics is taken into account the system gets destabilized. In fact for frequencies much faster than 
H and ciJgrad our system just reduces to the ghost condensate in Minkowski space, which features a 
Jeans-like instability [7j. In this regime we can estimate the rate of the gravity-induced instability 
in our system to be the same as that of the ghost condensate, 

2 

which is faster for larger M . In fact we get Wgrad "-i^ Jeans ~ -ff . In this case it is obvious that the best 
we can hope for is a compromise: the slower we make the gradient instability, the faster the Jeans 
instability becomes, and vice versa. As we discussed in section [2] we want both kinds of instability to 
be slower than Hubble. Therefore the parameters M and M must be chosen so that the combination 
MM'^/MI^ hes in the interval 

— < ^Jli^ < H (47) 

This is really an interval only if H is parametrically smaller than H^. For H ~ the choice is 
highly constrained: MM'^/M^,, ~ H. Finally, for H much larger than one of the two instabilities 
is unavoidable. 

So far we neglected a possible term 

AS^atter = - [ d^x ^ M'^SEUSN . (48) 
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Given the trasformation laws of 6N and 6Eij, in terms of the vr field this operator corresponds to a 
three-derivative term of the form — M^vr^V^vr, apart from a negligible correction to the vr^ term of 
Eq. (j38p . However upon integration by parts such a term can be rewritten as 

- M^T^lv^TT ^ M^^^^{di7rf ^ -HM^^id^Trf , (49) 

where again we are assuming that the time variation rate of the coefficient M is generically of order 
H. We therefore get a gradient energy term which, although it is suppressed by H, can in principle 
compete with that coming from the tadpoles, Eq. ([36]) . itself suppressed by H. In particular by 
choosing larger than roughly HMpjH we can make the gradient energy positive. But again we 
have to worry about the effect of mixing with gravity. In analogy with the above discussion such a 
mixing is unrelated with the background stress-energy tensor and thus with H, and gets enhanced 
when we take larger and larger. This makes the Jeans-like instability — if present — faster and 
faster. However there is a qualitative difference with respect to the previous case: here there are 
no terms, so that the mixing with gravity is relevant not only in the IR but at all scales. It is 
easy to check that the mixing terms become as important as the diagonal ones for > HMpy As 
we argued, the gradient instability is cured for > HMp^/H. Then for H <C one can choose 
M^/Mpy to lie in the same parametric range as before, Eq. (j47p : in such a range mixing with gravity 
is negligible and the system is stable at all wavelengths, with no need of higher-derivative terms in 
the UV and of the cosmological horizon in the IR. As we will see in the next section, systematically 
taking into account the full gravitational dynamics shows that this parametric range for is indeed 
the only one for which the system is stable. 

In conclusion — as long as we restrict ourselves to systems that spontaneously break time trans- 
lation invariance — we verify the generality of the tension found in section [2] between a violation of 
the null energy condition and the stability of the system. The system can be made completely stable 
only when H <^ H^. We are now going to confirm the results we got by a full calculation including 
gravitational effects. 

3.3 Full gravitational analysis 

In this section we want to write an explicit second order Lagrangian for scalar perturbations around 
an FRW background with H > including the mixing of the scalar with gravity. This will show 
both the instabilities we discussed and the constraints on the model parameters to achieve stability. 

Notice that the ADM variables N and iV* are not dynamical, i.e. the action does not contain 
their time derivative. They should be thought as Lagrange multipliers: their equations of motion are 
respectively the Hamiltonian and momentum constraints of General Relativity. This means that we 
can solve for and from their equations of motion and plug the result back into the action, to 
get the Lagrangian for the scalar mode we are interested in. The Lagrangian turns out to be more 
useful to study instabilities than directly looking at the linearized equations of motion. A ghost 
instability, where a system has the "wrong" sign of the energy, does not show up if we study the 
linearized equation of motions, but it will be manifest at the Lagrangian level. 

We choose unitary gauge, vr = 0. We have still to fix space diffeomorphisms and a convenient 
gauge is given by 

hij = a{tf [(1 + 2C)5ij + 7ij] , daij , 7ii = . (50) 



15 



The matrix 7 describes tensor modes and in the fohowing it will be neglected as we are only interested 
in the dynamics of the scalar perturbation In this gauge the scalar perturbation is given by an 
isotropic perturbations of the 3D metric at constant t. 

Let us now proceed to write an explicit Lagrangian for at second order. Postponing the 
discussion about the SN5E^i term, we start from the action we constructed in the previous section, 

S = I d?xdt ^ ^ + N-\EijE'^ - E\^)) - Ml, (^-^H{t) + 3H{tf + H{t)^ 

(51) 



^^SN^-^5E^^ 
2 2 



In the first line the first two terms reconstruct the full 4D Ricci scalar, as in Eq. (j24p . The other terms 
of the first line are the tadpole terms as in Eq. ()35p . while on the second line we have contributions 
which are quadratic in the perturbations. For simplicity we neglect quadratic terms of the form 
6Eij6E^^ , as we explicitly checked that they do not change the results, and we take M and M time 
independent. 

The variation of the action with respect to N gives the Hamiltonian constraint 

- N-'^{E,jE'^ - - i^-N-^H + 3H'^ + Hj + M*6N = , (52) 

while the variation with respect to gives the momentum constraint 

MliV, {N-\Eh - SjE'^k)) - M^ViE\ = . (53) 



We are interested in the Lagrangian at second order. Thus we have to solve these equations only at 
first order in the perturbation ^, as second order terms would multiply, once substituted back into 
the action, the unperturbed constraint equations dL/dN and dL/dN^, which vanish. Expressing 
the perturbations as = 1 + A'^i and A^* = diip we can solve the two constraint equations in terms 
of C 



{-ISMl^M^H'^ + 2M2jM4 - AM^^H) ( - AM^^H ■ V^C/a^ 



(54) 
(55) 



where we assumed M, M <C Mpi. 

The quadratic action for ^ is then obtained by substituting these solutions into the original action 
Eq. ()5ip . After some integration by parts, the result is given by 



S 



(fxdta-^{t) 



(56) 



We have three operators quadratic in C, with coefficients which are time dependent and given by 



A{t) 



2 M^i (m-^ - 9M^ H'^ -2 M|i H 
4M4j m + M'^M^ 



(57) 
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Bit) = £^i=T5 -24M^,M + M (M^ - 2Mi^H){M^M - 4.M^^H) 

+ 4ij2(M'^M|iM^ + 4M|ii7) - 8Ml.yM^HH] (58) 

^ 2M$,M^ 

4M^^H^ + M 

Using the ADM approach we managed to get a Lagrangian for the single relevant degree of freedom 
Q. It contains all the information about the dynamics of the system, in particular all the possible 
instabilities we are interested in. Notice that we have no mass term for ^, but only terms containing 
derivatives. The reason is quite general and it extends beyond the quadratic order [16]. From 
the definition Eq. (j50p it should be clear that if we can neglect spatial gradients of Ci i-e. for a 
sufficiently long wavelength, this variable cannot evolve in time, as it is equivalent to a constant 
isotropic rescaling of the spatial coordinates. A mass term would not allow this and it is therefore 
forbidden. Notice that the action for Q allows, disregarding gradient terms, to study the homogeneous 
perturbations of the background and in particular to check if a background solution is an attractor. 
From the time kinetic term of C, in its action Eq. ()56p we immediately get the solution 

In this way, it is easy to check that slow-roll inflation and the ghost-condensate are attractors. 

Let us now proceed to explicitly recover from this Lagrangian all the qualitative results about 
instabilities of the previous sections. As discussed, we are interested in time scales much smaller than 
the Hubble time. This simplifies the algebra; for example in deriving the equations of motion from 
the action above we can neglect the time variation of the coefficients A, B and C . An additional 
simplification comes from the fact that for the validity of the effective field theory description we 
must have H <^ M and H <^ M. We thus reduce to a rather simple equation of motion 

where k is the physical momentum. As we discussed in the previous sections, in Eq. (j57p we see 
that to avoid ghost-like instabilities we have to take > 2M^^H. We can focus on the regime 
3> 2MpyH, because if the two terms are comparable the instabilities are clearly worse. The 
resulting dispersion relation captures all the qualitative features we discusseci§ 

2 -{m'^M'^ + AM^,H)k'^ + 2M^,M^k^ , , 

" = 2MPf^ ■ ^^'^ 

The two terms proportional to /c^ describe the two instabilities we have discussed: the first, which 
is independent of i^T, corresponds to the Jeans instability already present in Minkowski space, while 
the second leads to the gradient instability, weighted by H. It is straightforward to find the unstable 
momentum k with the fastest (imaginary) frequency; this defines the instability rate. We get 

nUM^ ■ Ml \ , , 



^In the limit we are discussing one can check that the inclusion of the SEijSE''^ term just amounts to a 
redefinition of the constant M. 
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This formula unambiguously displays the complementarity of the two physically distinct instabilities 
we discussed; neither for MM^/Mpj — nor for MM^/Mp^ ^ oo do we end up with a stable system. 
One can try to minimize Wmax by a proper choice of MM^/Mp^. The rate of instability is minimized 
for 

MM'^/M^y = 2H^/^ Lo^^ = . (64) 



The full analysis therefore confirms what obtained in the previous sections. The instabilities can be 
pushed out of the horizon by a proper choice of the model parameters only if 



(65) 



Let us now analyze the effect of the SNSE^i operator. For simplicity we set to zero the SE^'i'^ 
term, and consider the following action 



S 



J 



+ 



(fxdt 



-9 



Ml, 



+ N-\EijE'^ - E\^)) - Ml, (^j^H{t) + ^H{tf + H{t) 



SN^ - M^5N6E' 



(66) 



The calculation follows exactly the same lines as before and we end up with an action for C, of the 
form 

a 



A{t)C,^ + B{t){-0' 



S = j d^xdta^{t) 
where the the time-dependent coefficients are 

Ml,{-\2Ml,HM^ + 2MUf|i - AM^^H) 



(67) 



Bit) = 



(-M3 + 2Ml,Hf 

Ml,{M^' - 2Ml,HM^ + 4M|iiJ) 
(-Af3 + 2Ml,Hf 



(68) 



(69) 



Notice that, as expected, we have no terms, which in the previous case were coming from the 
operator SE^i"^. In the limit of time scales much smaller than we get the dispersion relation 



a;2 = 



-M^ + 2Ml,HM^ - AM^^H 



'PI 



(70) 



Clearly the system is stable only if the coefficient of the k'^ term is positive. We can obtain this with 
a proper choice of M, only if 

■ (71) 



H < 
~ 4 

In analogy with the previous case, there is a parametric window for M only if ij ^ H"^ 



(72) 



In conclusion, the explicit calculation confirms what we argued in the previous section. Both 
with the SE'^iSE^i and with the SNSE^i terms we can have stability only if ij < H"^. In a generic 
effective action both quadratic terms will be present at the same time but at this point it is clear 
that this cannot change the qualitative picture. 
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4 Applications 



In this section we give some illustrative examples of non-standard cosmological histories that are 
now allowed as the NEC can be violated without introducing pathologies. 

4.1 Today's Acceleration 

After the surprising discovery that the universe is presently accelerating (in the sense that a > 0), 
it is natural to consider the even more exotic possibility that the Hubble rate is growing with time 
{H = a/a — (d/a)^ > 0). It has become common to parameterize the present expansion in terms 
of the equation of state parameter w = p/p, where p and p are respectively the present pressure 
and energy density of the universe. The violation of the null energy condition is equivalent in this 
language to the inequality w < —1. Present data require w > —1.2 [Il[2], leaving room for observing 
w < -1. 

The simple toy model in section [2] gives an example of a model with w < —1 with no fine tuning 
or instability. More general models can easily be written using the formalism of section [3l As already 
discussed, to avoid instabilities the model parameters must satisfy the inequalities 



more interesting to consider the case where w is near the lower experimental limit |w + 1| ~ 0.2. 
In this case, the instability is not parametrically suppressed, and we expect that the new degree of 
freedom has an unstable mode. We stress again that this is a long-wavelength instability like the 
familiar Jeans instability for matter, and not a breakdown of the effective theory. The presence of 
unstable modes in dark energy in models with -|- 1| ~ 1 is very general, and gives a possible new 
observable handle on these models. We leave the details for future work. 

We want to stress that this way of explaining a possible observation ofw < —1 relies on a genuine 
violation of the null energy condition: the Einstein frame metric has H > 0. This is completely 
different from other "phenomenological" approaches that mimic a violation of the null energy con- 
dition. For example one way to get < — 1 is to couple matter to a conformally rescaled metric 
f{(p)gfj,u, where gfj^i, is the metric in Einstein frame and (p is an evolving scalar field. Experiments 
will be sensitive to the way this new metric evolves, so that an observation of < —1 would not 
imply a violation of the null energy condition. Notice however that this kind of models is severely 
constrained by fifth force experiments, see for example Refs. [17^ [6]. Other approaches [181119] mimic 
the presence of a dark energy component with u; < — 1 while keeping H < 0: the total equation of 
state of dark matter and dark energy has u; > — 1. 

4.2 Starting the Universe 

If the NEC is satisfied, a very general property of an expanding universe is that it always evolves 
from a state with high energy density towards a state with a lower one. This implies that any 
effective field theory description will eventually break down if we go sufficiently backwards in the 




(73) 



or analogous ones in the presence of a 6N6E^i term. If we demand that instabilities are para- 
metrically suppressed, we must have H <C H^, and therefore -|- 1| <C 1. It is phenomenologically 
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past. Now that we have shown that it is possible to violate the NEC without introducing patologies, 
we can consider the possibility that the universe "starts" from a very low energy state, gains energy 
in the expansion, and eventually reaches an high energy state from which the standard cosmological 
evolution begins. In particular, the universe could approach flatness in the asymptotic past, H ^ 
for t — > —oo. 

Notice that in any such model, independently of the details of the evolution, the present causal 
horizon is always infinite, since the universe expanded for an infinite proper time. This gives a 
potential alternative approach to standard inflation for solving the homogeneity problem of the 
universe. 

It is straightforward to implement this in the framework of the general effective action presented 
in section [3l specifying for example in Eq. (|4ip a suitable function H{t). For instance, if we take 
the scaling solution H{t) ~ a/\t\ for t — oo with a 3> 1, we approach flatness in the past, while 
keeping the required parametric separation 

„2 



^-^«^-^'- (74) 



The instabilities are under control if 
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H^ Mit)Mit) 

Ml, - ^ ^ 

Since H{t) for t — > — oo, we see that the model parameters must change with time in order 
to keep these conditions satisfied. As we discussed, since the system is not time-translationally 
invariant, it is generic that the coefficients of the effective Lagrangian M, M, M, . . . explicitly depend 
on time with a time scale of order H. There is therefore no difficulty in making the system stable 
at all times. 

It is easy to give an explicit realization of this scenario in a scalar field model similar to that 
of section [2j Besides an appropriate choice of the potential, an explicit cp dependence must also be 
present in the derivative terms. At the level of the vr action, this translates into a time dependence 
of the parameters, so that with an appropriate choice, the condition Eq. ()75p can be satisfied. Notice 
that the explicit (p dependence is compatible with the fact that the shift symmetry on (j) is already 
broken by the potential. To give an example, we start from a Lagrangian of the form Eq. ([3]), which 
admits the solution (j) = t. We choose the minimum of P to be at zero, e.g. 

P{X) = i (X - 1)2 , (76) 

so that in the absence of a potential term the solution for the metric is Minkowski space. Since we 
want H to vary slowly with time, the stress-energy tensor must be dominated by the potential. To 
reproduce H = a/\t\ we need 

V -/i . (77) 



As in section [21 this will slightly displace (p from the minimum of P, but this perturbation is small at 
early times and the corresponding correction to the stress-energy tensor is negligible. This solution is 
approaching a "big rip" singularity as t — > 0, so at some negative t it will break down. In particular, 
at t ~ —a we would have ttq ~ 1, and the unperturbed solution cannot be trusted any 
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more. This corresponds to a maximum Hubble rate of order 



max 



~ a 



/^mVMpi. Before this 



happens, we can match the solution to a standard FRW phase. 

To complete the construction of the model, we must introduce higher derivative terms. To satisfy 
the condition Eq. (j75p at all times, we need M to depend on time, which is accomplished by making 
the 4-derivative terms depend on (/> as |§ 



One can check that the contribution of this term to the stress-energy tensor is negligible at sufficiently 
early times, up to Hubble rates of order M^/Mpi. The (/)-dependence of the potential and of the 
(□(/))^, besides making the mass scales time-dependent, introduces small changes in the dynamics 
of the perturbations. For example, masses and additional (9vr)^ terms will be generated, but they 
will be negligible. This can be simply understood from the fact that for large a the rise of H is very 
slow and the system approaches the ghost condensate in de Sitter space. 

In summary, putting together the Lagrangian terms Eqs. (j76H78p we have an explicit model for a 
universe which starts expanding from a zero-curvature state without developing any instability. The 
point is not to present a particular compelling cosmological model, but to stress that it is possible 
to build models where the present universe evolved from a state which is asymptotically flat in the 
past. 

4.3 The Eternally Expanding Cyclic Universe 

Given the two scenarios above, one is tempted to put them together to build an eternally expanding 
universe with a cyclic evolution. We can imagine, as in section WA\ that the present expansion has 
H > 0. This means that in the future the universe will be inflating with a larger and larger energy 
density. At a certain point this energy can be converted to matter and radiation as in a conventional 
inflation model to start a new FRW- like evolution |20l [TO] . Then one is lead to connect the present 
acceleration with the inflationary phase which occurred in our past, to build a periodic universe that 
goes through this cosmological history many times. 

An illustrative example is based on the model of section [2j At the present epoch, the field rolls 
up the potential giving a super-accelerating phase with growing energy. The reheating to a new 
FRW phase can be achieved by a sudden drop of the potential. If we imagine that the field (p is 
periodic (see Fig. [2]), we obtain a periodic cosmological evolution in which the present acceleration is 
the beginning of a super-inflationary phase identical to the one that is responsible for the structure 
formation in our universe. 

As already discussed in section [2l the linear potential is technically natural because the shift 
symmetry is unbroken in the absence of gravity. The sharp drop of the potential responsible for 
reheating breaks the shift symmetry even in the absence of gravity, and will therefore induce more 
general terms that break the shift symmetry through radiative corrections. However, the breaking 
of the shift symmetry is localized in (p, so it is technically natural for the potential to be linear a 
distance A(/) > M away from the sharp drops. 

^In this way we introduce a large hierarchy between M and M which may be difficult to justify at the 
level of effective field theory. Alternatively both scales could move together remaining in the allowed range 
Eq. ([75)) . In this toy model, we stick to the simplest case in which only M is time dependent. 




(78) 
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The model can be made stable for the entire cosmological history. If at the present epoch Eq. (j73p 
is satisfied, we know that we are safe from instabilities while the scalar rolls up the potential. Looking 
backwards in time from the present epoch, we enter into a matter dominated and then radiation 
dominated phase (see Fig. [2|). As is larger then nowadays, the Jeans instability, whose rate is 
independent of H, is under control. Also the gradient instability is not an issue since its rate is 
proportional to V / H and it thus becomes smaller and smaller in the past. 




Figure 2: Schematic representation of the periodic model presented in section [43l 



If we take this model seriously, it gives a striking relationship between the present equation of 
state w and the tilt of the inflationary spectrum n^: they both depend on the constant slope of the 
potential V. 

Following Ref. [201 HH] the spectrum of density perturbations goes as 

n~(§)' P!.) 

where Hj is the Hubble rate during the inflationary phase. The slow variation of Hj gives a tilt: 

5 V 

n, - 1 ~ ^ — ^ . (80) 

On the other hand, the current acceleration of the universe would be characterized by an equation 
of state w equal to 

v 

l + w = ^ — T , (81) 

where Hq is the current Hubble rate. We thus get the relation 

n.-l = -y [j^J (l + w). (82) 

Given the constraints on w, the tilt of the spectrum is predicted to be extremely small and practically 
impossible to measure. Moreover, to be able to test this relationship one should determine Hj. The 
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production of gravitational waves does not help, since the energy scale M is constrained to be too 
low. In principle M could be determined from gravitational experiments or through the direct 
coupling of this scalar sector to the standard model, so that Hj could be fixed through Eq. ([79]) . If 
the detection of a negative tilt by the WMAP experiment [1] is confirmed this specific model will 
be ruled out. 

If we allow the scale M to be a function of time we do not have constraints coming from the 
present day modification of gravity. In this case the amplitude of gravitational waves could be 
sufficiently large to be observed. Notice that these gravitational waves would have a blue spectrum, 
a striking signature of violation of the NEC. 

In conclusion we have shown how the violation of the null energy condition opens up the possi- 
bility of relating the present and past evolution of our universe and of building periodic cosmological 
histories. 



4.4 The Big Bounce 

The null energy condition implies that the cosmological evolution in the absence of spatial curvature 
cannot bounce from a contracting to an expanding phase, since we need -ff > to make H change 
sign. We expect that bouncing cosmologies can be realized in our scenarios. However, as -?/ = at 
the bounce, one could wonder whether the instabilities can be kept under control. Obviously the 
real issue is whether the instabilities have sufficient time to grow. We will see that it is possible to 
control instabilities by making the bounce sufficiently fast. 

Let us take the following simple evolution of the scale factor around the bounce (see figure [3]) 

ait) = 1 + ^ , (83) 

so that for t <^T 

Hit) ± . (84) 

A system that realizes this bounce is described for example by the vr action Eq. (j4ip . with Hit) 
as above. The stability condition Eq. (j47p shows the potential problem. The Hubble rate, which is 
crucial in stabilizing the Jeans and gradient instabilities, goes to zero at the bounce. 

Let us see how the problem can be solved. For simplicity we assume that the parameters M 
and M are time-independent. Notice that during the bounce H is constant, H ~ As we have 

shown in section O for a given H the instability rate is minimized for MM'^ /Mpy = 2ij^/^, which 
gives 

C^inst = H'/^ = ^ . (85) 

This choice makes the rate of gradient and Jeans instability the same; see Eqs. (|63H64|) . 

In this case the only time scale during the bounce is T. In particular at t ~ itT the Hubble rate is 
iJ(ibT) ~ ibl/T. Therefore the period for which the instability is faster than H only lasts a time of 
order T, which is one instability time. That is, the instability does not have time to develop. There 
is, however, another form of instability one has to worry about during the contracting phase. In a 
contracting universe anisotropic contributions to the Friedmann equation blue-shift like a~^. This 
means that in order for a contracting solution not to be destroyed by anisotropic metric perturbations 
the contraction must satisfy H < —3H^ [21]. Of course when the bounce phase sets in at t ~ —T 
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Figure 3: Schematic dependence of the scale factor a{t) and of the Hubble rate H{t) in the bouncing 
model of section assuming that the bounce is matched to standard H < phases both for t < —T 
and for t > T. 



this condition is violated, but since the scale factor only varies by a factor of order unity during the 
bounce, the anisotropics will only grow by an order one factor, which obviously is not a problem. 

Outside the |t| < T window, one can match the bounce to a contracting phase for t < —T and 
to an expanding one for t > T. The expanding phase can be a standard decelerated FRW evolution, 
similarly to what happens in ekpyrotic/cyclic scenarios, or a super-accelerating phase with H > 0, 
which can make H as large as we like. For negative t one is forced to consider a H < phase to 
avoid growing anisotropic perturbations. Notice that as H decrease going to negative t (and also for 
positve t ii H < 0) the Jeans instability is active. This can be cured as we did in section 14. 2^ e.g. by 

2 

postulating that M has a suitable time-dependence that makes the Jeans instability rate slower 
than H. This is certainly consistent at the level of the effective field theory for the fluctuations 
we developed in the previous section. In the language of the scalar cj), this would correspond to a 
particular (p dependence of the higher derivative part of the Lagrangian. Finding the explicit form 
of such a Lagrangian is out of the scope of the present paper. 

We conclude that it is possible to realize a FRW solution in which the universe bounces from 
a contracting phase to an expanding one, without developing catastrophic instabilities. To our 
knowledge, this is the flrst case in which a bouncing solution is realized in a situation where the 
importance of all the higher derivative terms is controlled in an effective field theory expansion at 
low energies. 

In this controllable setup it is easy to follow the evolution of cosmological perturbations from the 
contracting to the expanding phase. This is crucial for the experimental viability of ekpyrotic/cyclic 
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scenarios [9j. Since our bounce satisfies the general hypotheses of Ref. [22], it unfortunately leads 
to a non scale- invariant prediction for density perturbations. 

5 Conclusions and Outlook 

In this paper we have shown that at the level of effective field theory there is no obstruction in 
constructing systems that violate the null energy condition (NEC). In particular, the instabilities 
that plague a very broad class of NEC-violating systems [6] can be completely avoided. We use this 
to construct cosmological models with H > that do not have any instabilities. 

These results depend crucially on the interplay between two potential types of instabilities. The 
first are gradient instabilities that occur very generally whenever the NEC is violated. These can 
be avoided due to higher-derivative terms in the effective theory. However, increasing the higher- 
derivative terms increases the effect of another potential instability, the Jeans instability that arises 
from the attractive nature of gravity. The Jeans instability can be avoided if the time scale for the 
instability is longer than the Hubble time. Demanding that both of these instabilities are absent 
requires 

H<,H\ (86) 

This means that we cannot get an arbitrarily large violation of the NEC. These results are derived 
both in a specific model (a deformation of ghost condensation [7] deformed by the addition of a 
potential |20[llUj) and in a very general effective field theory analysis of adiabatic scalar fluctuations 
about an arbitrary FRW background. 

We emphasize that if Eq. (|86p is violated the effective field theory does not break down, but it 
contains long wavelength exponentially growing modes. These might be interesting for cosmology: 
for example, if the present acceleration of the universe has w = p/ p < —1, then it becomes possible 
that they might have shown up around our present epoch. 

We constructed a number of explicit examples to show how this can lead to interesting alternative 
cosmological scenarios. We constructed a simple explicit model that gives w < —1 today. We also 
presented a model in which the universe starts from Minkowski space in the distant past, giving a 
possible alternative to standard inflation for solving the horizon problem. Putting these together, 
we consider an eternally expanding cyclic model in which the present accelerating phase is the 
beginning of inflation in the next cycle. This opens up the intriguing possibility that measurements 
today can give information about the inflationary phase that gave rise to the structure we observe 
today. Finally, we constructed a model which has a smooth "bounce" connecting a contracting phase 
to the present expanding phase. Although bouncing cosmologies have been previously considered in 
the literature, we believe ours is the first example in which the bounce is under theoretical control. 
These illustrate some of the phenomenologically interesting possibilities opened up by violation of 
the NEC. We leave detailed analysis of these ideas to future work. 

We close with some theoretical questions. It would be interesting to understand whether the 
requirement H <C to have a parametric suppression of instabilities is completely generic or not. 
Models involving more degrees of freedom could be completely stable even when the NEC is strongly 
violated. 

While the results of this paper are under control in a systematic effective field theory expansion, 
it is natural to ask whether this particular effective field theory can be embedded in a UV complete 
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theory such as renormaUzable quantum field theory or string theory. In particular, the somewhat 
exotic features of the ghost condensate (e.g. Lorentz breaking and complicated nonlinear dynamics 
[7J) may make one suspicious that the ghost condensate lies in the "swampland" of effective field 
theories without UV completions [23^ I24j. In fact, in Ref. [25] it was shown that a large class of 
consistent effective field theories (including some interesting modifications of gravity) have no Lorentz 
invariant UV completion, essentially because they allow configurations in which signals travel faster 
than light [25]. The present model does not suffer from this problem, since the excitations of the 
ghost condensate have a dispersion relation uJ^ ~ k^/M"^, and therefore travel much slower than 
light. However, it is true that there is currently no known UV completion of the model; for an 
interesting recent attempt, see Ref. |26j . 

Despite these open questions about the full consistency of the model at all energy scales, we 
believe that the phenomenological possibilities opened up by violation of the NEC are well worth 
further exploration. 
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Note added 

In a recent paper [27j a NEC-violating model was presented. It involves a vector field as well as 
a scalar. Like ours, the model has no UV — i.e., arbitrarily fast — instabilities, but only exponential 
instabilities below a critical momentum kc- However, not even in the 'slow roll' regime H <^ can 
the instability rate ujc be made slower than the Hubble rate. What prevents the system from being 
disrupted is then that fluctuations exit the horizon, so that instabilities only have a finite amount 
of time to develop before being damped by Hubble friction. On the other hand, in our case when 
H <^ the system can be made completely stable. 
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